Electron-like and photon-like excitations in an ultracold Bose-Fermi atom mixture 
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We show that the electron-like and photon-like excitations may exist in a three-dimensional Bose- 
Fermi Hubbard model describing ultracold Bose-Fermi atom mixtures in optical lattices. In a 
Mott insulating phase of the Bose atoms, these excitations are stabilized by an induced repulsive 
interaction between 'electrons' if the Fermi atoms are nearly half filling. We suggest to create 
'external electric field' so that the electron-like excitation can be observed by measuring the linear 
density-density response of the 'electron' gas to the 'external field' in a time-of-fiight experiment of 
the mixture. The Fermi surface of the 'electron' gas may also be expected to be observed in the 
time-of-flight. 
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Ultracold atoms in optical lattices have offered a highly 
tunable platform to study various physical phenomena 
which may not be definitely clarified in condensed mat- 
ter systems p]. On the other hand, new systems which 
may not be realized in condensed matter content are pre- 
sented, for example, mixtures of Bose-Fermi atoms as 
constitution particles. Experimentally, the Bose-Fermi 
atom mixtures in optical lattices have been realized for 
87 Rb- 40 K 0, and 23 Na- 6 Li 0. 

Microscopically, these mixtures may be described by 
a Bose-Fermi Hubbard model 4]. The constitution par- 
ticles in this model are spinless boson(ai) and fermion 
(/i) with i the lattice site index. In this Letter, we con- 
sider three-dimensional cubic lattices. We will show that, 
in high temperatures, there are only excitations of these 
constitution particles. We call this a confinement phase. 
The system undergoes a phase transition, in certain crit- 
ical temperature, to a uniform mean field (UMF) state 
in which electron-like and photon-like excitations emerge 
if the fermion occupation is nearly half filling . At the 
exact half filling, this UMF state turns to a long range 
ordered state, the checkerboard crystal of 'electrons' |g. 
We show that this UMF state can only be stable if the 
bosons are in a Mott insulator (MI) ground state. 

We suggest an experiment to create an 'external field' 
by changing the depth of the fcrmion's optical potential 
7]. The response function of the mixture to the external 
field may be measured by the density distribution image 
in a time-of-flight of the mixture cloud. The behavior 
of the response functions may be used to identify the 
electron-like excitation. We expect the 'electron' Fermi 
surface can be observed by the time-of-flight experiment, 
which has been used to observe the Fermi surface of the 
pure cold Fermi atoms 0- 

The Bose-Fermi Hubbard Hamiltonian we are inter- 
ested in reads 

H = -J2( t Ba\a J +t F flf J )-J2^Bnt + li F n{) 

(ij) ' 
U BB 



where the lattice spacing A/2 is set to unit. ( We also 
set h = c = 1.) = a\a,i and n{ = ts and tp 

are the hopping amplitudes of the boson and fermion be- 
tween a pair of nearest neighbor sites (ij). \ib and A*f 
are chemical potentials. And Ubb and Ubf are the on- 
site interactions between bosons, and between boson and 
fermion. In this work, we use Ubb > Ubf > although 
this is not necessary in general. The microscopic calcu- 
lations of these model parameters in terms of the cold 
atom mixture have been established, e.g, in Ref. 0. 

To deduce the low energy theory in strong couplings, 
we will use the slave particle technique, which has been 
applied to the cold boson system j^, [nj uM ■ ^ n the slave 
particle language, the Hamiltonian reads H = H% + H4 
where H2 and H4 are the two-operator and four-operator 
terms, respectively Namely, 
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where = c' a i c a ^ and = h' aji h a ,i- We explain 
briefly the derivation of this slave particle Hamiltonian. 
The state configurations at an arbitrary given site con- 
sists of {\a, s) I a = 0, 1, 2, ... ; s = 0, 1} where a and s 
are the boson and fermion occupations, respectively. The 
Bose and Fermi creation operators can be expressed as 



it = 



J2 a V^ + T[|a + 1, 0) (a, 0| + \a + 1, 1) (a, 1|], /t 



J2 a \a, l)(a,0|. The mapping to the slave particle reads 



fermion (CF) |f| and h, 



We call c a the composite 
the slave boson. The nor- 
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malized condition J2 a (\a, 0)(a, 0| + \a, l){a, 1|) = 1 im- 
plies a constraint ^2 a (n" + n^) — 1 at each site. The 
slave particles arise a U(l) gauge symmetry c Qj i — > 
e zipi c a ,i,h a ^ — ► e lVi h a ^. The global t/(l) symmetry 



* v h a ^ reflects the particle num- 



ber conservation. Since the slave particle technique es- 
sentially works in the strong coupling region, we focus on 

u BB /(m B ) > i. 

There are two types of four slave particle terms in 
if 4 , the t B -terms and ip-terms. We first neglect the 
is-terms in the mean field level of the CF. To decou- 
ple the i^-terms , we introduce Hubbard-Stratonovich 
fields ffep f . and Xa'p y ■ The partition function is given 
by Z = J DxDfjDcDcDhDhDA e~ S 'is where the effec- 
tive action reads 

f 1/T f 

Seff K> Vi c, h, A] = J dT^H 2 + 22 iA i 
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+ E t F $p a Ji(Xa0,ij ~ c a.iCf3,j) ~ Xa/3,ijXpa,. 

+vtp,ij(Xp a ,ji ~ hpjha,i)] + t B tern 
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where A, is a Lagrange multiplier for the constraint 
EaKi + = 1- Rewriting ffc^ = rf^^e^ , 

Xa/3,ij = Xap,ij eiAii and A 2 = A + Ao,i, Aoi and ,4# are 
£7(1 ) gauge field corresponding to the gauge symmetry 
[l2|. Before going to a mean field state, we first require 
the mixture is stable against the Bose-Fermi phase sep- 
aration. For example, it was known that the mixture 

,A7T__ ^ 7T2 



is stable in the MI phase if intp sin(7TO'' )U BB > U BF 

0. Near the half filling of Fermi atoms, this condition 
■ ■ - - — • h 

A (which is a saddle point value 
of Aj). This is corresponding to a mean field approxima- 
tion. The effective mean field action is given by 



is easy to be satisfied. Fixing a gauge, f]ap,ij " "'' 
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frequencies. Near the critical point where Xdp % 
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ons can expand -Djjy by x a 
mean field equations are X a p,ij = tFXp a 



For a j3, the 
\- and 
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tFXp^^jt with n% = [e-^B - l]-\n« = 

- The solutions of these mean field equa- 

tions can only exist in the weak coupling limit, i.e., 
ftps <C tp: — tp. Therefore, in the parameters 



we are considering, \ a 



for a ^ (3. For a = fi, 
the mean field equations are T^x a a,ij = T FX aa ,ji an d 
T?xL,ij = TgxZaji with T^ B = t F nl h {l T <J. The 
critical temperatures are then given by T" = ^/TpT B . 
Below T", the minimized free energy including variables 
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is F cx -ifEJ^) -gs- — ■ - -^r 
S%, S" are the fractions of non-zero terms corresponding 
to the second order and four order of \x aa \ 13]. The 
optimal \ X c aa | 2 ar tt S 2 a /5 4 a 

We consider an integer boson filling n a =1 for J7b.b 3> 
1. Since n-^ 1 and n"^ 1 in this region are very small, 
there are no solutions of the mean field equations with 
jtq^i > q f or ^- ne cr j^i ca i temperature equations. Thus, 
all slave bosons and CFs with a/1 are confined. Only 
T 1 > can be found. Below T 1 , the CF ci and the slave 
boson h\ are deconfined in the mean field sense. In the 
inset of Fig. 1, we plot for a set of given parameters. 
The curve T p (U B b) is corresponding to /j, F (T p ) = 0, i.e., 
the effective chemical potential of ci vanishes. 

We now go to concrete mean field solutions and fo- 
cus on the near half filling of the fermions. Because we 
work in a three-dimensional lattice, the flux quanta pass- 
ing a cubic cell are zero. Thus, there is no a flux mean 
field state. (The flux mean field state may exist in a 
two-dimensional Bose-Fermi Hubbard model.) Two pos- 
sible solutions are the dimer and uniform phases. In a 
cubic lattice, each lattice site has six nearest neighbor 
sites. A pair of slave boson and CF located at the near- 
est neighbor sites may form a bond. The dimer phase 
means for any given site, only one bond ended at the 
site is endowed with a non-zero Xaa an( l other five bond 
carry Xaa = 0- l n l ne uniform phase, each bond is en- 
dowed with the same real values of Xa h if a CF (i.e., 
a lfi\Q) = c iJ vac )) i s surrounded by six slave bosons 
(e.g., at|0) = h\-\va,c)). In the fermion half filling, this is 
a CF checkerboard crystal state |£g. Slightly away from 
the half filling, this is a UMF state with hi-ci bond. In 
the uniform phase, the dispersions of the CFs and slave 



( D h )ij == (°r - Mb)%<W - *F T,fX C a p,ji 6 j,i+r- with 

Ms = RbQ ~ ^y s -a(a — 1) — iA. (r is the unit vector of 
the lattice.) 

The mean field equations are Xapij = ( c a,i c P,j) = 

T T, n D c!ij(Pn),XaP,ij = ( h *,i h P,j) = T Y,n D h%(^) 

where D a c %{ Pn ) = Jdre^D^r) and D a h %(co n ) = 
J dTe luJ " T D^At) with p n and uj n the Fermi and Bosc 



At the fermion half filling, the mean field free energy 
favors for the dimer state. However, as we shall see soon, 
in the Mott insulator phase of the bosons, the boson hop- 
ping term (i^-term) we have neglected at the mean field 
level will contribute a nearest neighbor repulsion between 
CFs or slave bosons HUbb > U B f ■ This repulsion poten- 
tial will raise the energy of dimer phase while the uniform 
phase is not affected. Thus, for our purpose, we focus on 
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the uniform phase below. 

In the mean field approximation, we neglect the bo- 
son hopping term and the gauge fluctuations Aoi, Aij 
which must be considered if the mean field state is 
stable. We first deal with the is-term. Introduce a 
Hubbard- Stratonovich field <E>i = J2 a V a + l[c„ +1 iC a ,i + 

ha + i ih a ,i] to decouple is-term. $i may be thought as 
the order parameter field of the Bose condensation. The 
phase diagram of the boson may be determined by mini- 
mizing the Landau free energy associated with the order 
parameter (0|$i|0). The vanishing of the coefficient of 
the second order term in the free energy gives the phase 
boundary 0, . It has to point out that here the fluctu- 
ation Aoi has been neglected and a cut-off to the type of 
the slave bosons has to be introduced. However, the ex- 
perience to work out the pure Bose phase diagram showed 
that these approximations could be acceptable |9l llClllH . 

As expected, the phase diagram of the constitution 
boson consists of the Bose superfluid (BSF), the normal 
liquid and the MI, in which the MI phase only exists in 
the zero temperature and an integer boson filling factor . 
In the inset of Fig. 1 , we show the boson phase diagram in 
the same parameters as those in the CF mean field phase 
diagram. In the Bose condensate, the boson number in 
each site is totally uncertainty. This means the vanishing 
bond number or S2 ~ 1S4 ~ 0. Thus, the mean field 
state is not stable in the BSF. In the MI phase of the 
bosons, on the other hand, the boson number in each site 
is exactly one for n a = 1. Thus, the mean field states may 
be stable. Furthermore, the ^-term induces a nearest 
neighbor interaction between CFs or slave bosons. This 
may be seen by taking the ts term as a perturbation if 
Ubb and Ubb — Ubf are much larger than ts- To the 
second order of the perturbation, the is-term contributes 
an effective repulsive interaction between the CFs or slave 
bosons in the nearest neighbor sites 



'■J 



with J = 



16t|[fj, F 
Ubb(U% r -U- 



const = Jn^Tife,- + const', (6) 

m 

For the dimer state, this re- 



pulsive potential contributes an energy J/2 to a pair of 
adjacent bonds. For the UMF state, if the fermion is in 
half filling, the checkerboard distribution of the fermions 
(then CFs) makes J no contribution to the energy. For 
a doping 6, i. e., slightly away from the half filling, the 
energy raises a small amount of the order J 5. Thus, the 
UMF state minimizes this induced nearest neighbor re- 
pulsion. 

We have shown the mean field states are not stable 
in the BSF. Therefore, we focus on the MI phase of the 
bosons below and examine the gauge fluctuations. The 
zcroth component Aio restores the exact constraint of 
one type of slave particles per site while Aij restores the 
gauge invariance of the effective action These fluc- 
tuations may destabilize the mean field state. To see the 
stability of the UMF state against the gauge fluctuation, 
one should integrate out h\ and c%. In the long wave 



length limit (k — > 0), integrating out hi first and keeping 
the Gaussian fluctuations of the gauge field, an effective 
action in continuum limit is given by 



S[c 1 ,A fl ]=TY^ J d 3 k 



c\{ip n - Mc + ie A )ci(7) 



2m, 



■ci(k a + ie Q A a ) 2 a + -A^AuTl^k^n) 



where eo = 



A^ = A^/eo, m c ~ l/(t F xiS) is the 

effective mass of CF. The coupling constant 2ne 2 l = J/2 
is a small quantity means that the Gaussian approxi- 
mation to the gauge field is reasonable. n^„ is the re- 
sponse function by integrating over h\. The density- 
density response function is given by [n^w, k)]^ 1 = 

[H^'Vfc)]; 1 + V(k) where V(k) = J£ a (l - k 2 /2) 
is the Fourier component of the interaction © and 

nf" 
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(oj, k) is the free boson response function. The 
current-current response function has a form Ilfj — 
(S i: j-'hkj/k 2 )Ilf + (hkj/k 2 )Il B ,. In the long wave length 
limit, the transverse part n|! = — pTI^g. The longitu- 
dinal part nf (k,uj) oc |(0|/ii|0)| 2 as k,u> — > 0. However, 
out of the BSF phase, h\ does not condense. Thus, A^ 
is a transverse field and the action (JTJ is very similar to 
electron coupled to a photon field. We, therefore, call 
Ci and A^ the electron-like and photon-like excitation, 
respectively. Since the hopping of the Bose atom in the 
optical lattice is short range, the repulsive interaction 
between 'electrons' is also short range. If it was possible 
to design the hopping of bosons ts,ij oc l/yfi — j|, the 
interaction between 'electrons' would be purely coulom- 
bic, i.e., V(k) oc 1/k 2 . The stability of the UMF state 
against to the gauge fluctuations may be seen after inte- 
grating out ci field. This leads to a CF response func- 
tion [H^K,fc)]-! = [H^K,^]- 1 + [U^ n ,k)}-\ 
The UMF is stable when n^(0, 0) > 0. In the long wave 
length limit, it is positive if 6 J > m 7r fcl where k F is Fermi 

momentum of the CF. Since l/m c oc tpS, this condition 
holds only when the fermion filling is slightly away from 
the half filling. 

We have shown the stability of the UMF state near the 
fermion half filling when the bosons are in the MI phase. 
We may figure out the phase diagram of the CF in Fig. 
1. The mean field phase transition temperature T]. is 
suppressed greatly to T G which is determined by the T c 
and T p . The dash line is the estimated crossover line 
from the 'electron' gas (the MI of bosons) to the Fermi 
liquid of the constitution Fermi atoms (the normal liquid 
of bosons as the incompressibility of the MI is gradually 
disappearing). 

The experimental implications of the UMF phase are 
discussed as follows. We consider the 'electron' response 
to an external 'electric' field, 'made' by a change of the 
lattice potential of the Fermi atoms. This technique has 
been used to study the excitation spectrum of atom su- 
perfluid in optical lattices 0- This disturbs the density 
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FIG. 1: The phase diagram of the CFs in the Ubb-T plane 
for n a = 1, Ubf = Ubb/2, t F /6t B = V&j, n f = 0.55. The 
solid curves are the critical temperature T^* after considering 
the gauge fluctuation and the shrink of the CF Fermi surface. 

of Fermi atoms with H' = — s £lii n \ { Pi- in a time-of- 
fllght experiment, the difference between disturbed and 
undisturbed fermion densities by external field is given 
by n/ (r) - n° f (r) = (^)|«J,(k = ^f)\ 2 Sn f (k = ^f) 
where n/(r) is the image after the time-of- flight and 
71/ (k) is the Fourier component of nf; t is the flying 
time, Wf is the Fourier component of the fermion Wan- 
nier function and to/ the Fermi atom mass. If T > Tj}*, 
the density response of the system is simply given by 
the free fermion one, and then <Jn/(k) oc n^ ''. When 
the bosons are in the MI phase (nf — 1), n{ = n 1 ^ im- 
plies 5n f (k) = Snl(k). Thus, for T < T}* , especially 
below the dash line in Fig. 1, 5n/(k) cx IToo(k, 0) with 

Hoo 1 = (noo 0) ) _1 + (n^^V 1 + V(k). Then, this differ- 
ence between the response functions may be measured in 



experiment. A better experimentally measurable quan- 
tity is the visibility V = ^[r°^)+n^(rl'°j where r max 
and r m i n are chosen such that the Wannier envelop is 
cancelled 01- The difference between the disturbed 
and undisturbed visibility may directly correspond to 
the response function because n/(r max ) + n/(r min ) w 
n? (r roax ) + n°Jr m i n ) in denominator. To directly see the 
image of the fermion cloud, one may use a magnetic field 
to separate the fermion cloud from boson cloud before 
recording the fermion image in the time-of-flight as split- 
ting components in a spinor Bose atom condensate [15J. 

The Fermi surface of pure cold fermions has been ob- 
served in a recent experiment Q by the time-of-flight 
experiment. When the mixture in the UMF state, as 
we have discussed in the previous paragraph, n./(r) = 
(^ i )|zD/(k)| 2 nJ(k = Therefore, it is expected that 

instead of the Fermi surface of the free Fermi atoms, one 
may observe the 'electron' Fermi surface of the 'electron' 
gas. Namely, in the image of the time-of-flight, most of 
Fermi atoms are inside of area with |r| < |k^|i/m/ but 
not |r| < |ki?|£/ m / ( |kp| is Fermi surface of the free 
Fermi atoms.). 

In conclusions, we showed that there may be electron- 
like and photon-like excitations in mixtures of Bose- 
Fermi atoms in optical lattices in which the bosons are 
in the MI phase (Ubb and Ubb ~ Ubf Ib) and the 
fermions are nearly half filling. (To avoid the demix- 
ing, Antp sn\(imf)UBB > U BF .) It was suggested that 
through the time-of-flight experiment, the electron- like 
response function and the 'electron' Fermi surface may 
be measured. An electron-like response function also im- 
plies a gauge field effect. However, an experiment how to 
directly observe a 'photon' was not designed yet. 
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